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Abstract 

We consider the time evolution of a particle subjected to both a uniform electrostatic field F 
and a one-dimensional delta-function potential well. We derive the propagator ICf{x, t\x' , 0) of this 
system, directly leading to the wavefunction ippix, t), in which its essential ingredient /Ci?(0, t|0, 0), 
accounting for the ionization-recombination in the bound-continuum transition, is exactly expressed 
in terms of the multiple hypergeometric functions F{zi, Z2, - ' ' ; -^n)- And then we obtain the ingre- 
dient /Cir(0, t|0, 0) in an appropriate approximation scheme, expressed in terms of the generalized 
hypergeometric functions pFq{z) being much more transparent to physically interpret and much 
more accessible in their numerical evaluation than the functions F(zi,Z2, ■ ■ ■ , Zn)- 
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I. INTRODUCTION 



The ionization of atoms in an external electric field is, as well-known, one of the oldest 
subjects in quantum physics. As a simple example of the ionization, hydrogen-like atoms 
in a uniform electrostatic field have been extensively considered (see, e.g., [HE])- Here the 
background potential caused by the field decreases without limit in one direction, and the 
electron initially in the bound state will eventually tunnel through the "barrier" created by 
the field, so leading to ionization of the atom. The tunneling rate for an ensemble of many 
independent electrons has been calculated based on the exponential decay law following from 
the statistical assumption that the tunneling rate is proportional to the number of available 
atoms. 

On the other hand, big experimental advances in the field of nano-scale physics have 
highly enlarged a need for a detailed understanding of the (time-dependent) tunneling pro- 
cess of individual elections subjected to an external field. In fact, the nano-scale devices have 
been designed and fabricated, examples of which are molecular switches j3] and resonant 
tunnel junctions etc. However, in exploring the time evolution of the electron tunneling 
process leading to ionization, we have a considerable mathematical difficulty that there are 
no exactly solvable models for a transition from a bound state to the continuum. Further, 
even obtaining the numerical solution, with high accuracy, to this problem cannot be con- 
sidered an easy task either, especially in the strong-field limit where a highly oscillatory 
behavior is found in the time evolution of the bound-continuum transition. 

The system under investigation in this paper is a particle subjected to an attractive one- 
dimensional delta-function potential, —VoS{x). We intend to study analytically the time 
evolution of the particle when a uniform electrostatic field F is applied. The delta-function 
potential well (not in an external field) has a single bound state. This aspect may rather 
simplify the analytical study of the ionization. However, no exact solution %l)p{x,t) to the 
time-dependent Schrodinger equation of this problem has been found in closed form (in 
terms of the actual calculability to any sufficient degree of precision) , even for such a simple 
form of external field. 

In fact, this model of the field-induced time-dependent ionization has already been studied 
by some people. It was first discussed by Geltman [5j. Later on, several different approaches 
to obtaining the time-dependent wavefunction ipp{x^t) have been carried out [6iillj. One 



2 



of the approaches mainly adopted so far is to turn the relevant time-dependent Schrodinger 
equation into the Lippmann-Schwinger integral equation [cf. Eqs. Q and ([T])]. This integral 
equation is, however, highly non-trivial to solve analytically and so has been focused mostly 
upon its numerical solvability. In jTl [T2] by Elberfeld and Kleber, interestingly enough, the 
time-dependent ionization probability has been investigated in the strong-field limit based 
on the numerical analysis of the integral equation. They also demonstrated numerically 
that the ionization probability obtained from the simple exponential decay law may be 
considered a fairly good approximation on the average to the exact result in the strong- 
field limit although this approximation, by construction, cannot account for the short-time 
ripples, observed in the exact one, which result from the ionization-recombination process 
in the bound-continuum transition [TJ [13]. 

The goal of this paper lies in a systematic derivation of the propagator }CF{x,t\x' ,0) of 
the system, directly leading to ipp{x,t) = f^^dx' ICp{x,t\x' ,0) ipoi^')- As our central find- 
ing, its essential constituent part /Cp'(0, t|0, 0), accounting for the ionization-recombination, 
is exactly expressed in terms of the multiple hypergeometric functions F{zi, Z2, - ■ ■ , and 
then in an appropriate approximation scheme in terms of the generalized hypergeometric 
functions pFq{z) being much more transparent to physically interpret and much more acces- 
sible in their numerical evaluation than the functions F{zi, Z2, - ■ ■ , Zn)- The general layout is 
the following. In Sec. [IT] we briefly review the well-known results of both the delta-function 
potential problem without an external field and the problem of a particle subjected to a 
uniform electrostatic field but not bound by the ultrathin potential well. In doing so, we 



also sophisticate the old results. In Sec. Ill, we derive an explicit expression of the propa- 



gator in terms of the multiple hypergeometric functions and then discuss its mathematical 



complexity. Sec. |IV| we obtain the propagator in approximation, expressed in terms of the 
generalized hypergeometric functions. Finally, we give the conclusion of this paper in Sec. 

El 

II. GENERAL FORMULATION 

The system under consideration is described by the Hamiltonian 

-Hp = ^ - Vo6{x) - xF{t), (1) 
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where an external field F{t) = F ■ Q{t), and Vq > is a strength of the (5-potential well. 
For the field-free case (F = 0), it is well-known that the Hamiltonian T-Lq has a single bound 
state [H] 

Mx) = 75 e-^ 1^1 (2) 

with its eigen energy Eb = —h?B'^ /2m where B = mVo/h'^. All eigenstates and eigenvalues of 
Tio as well as the completeness of the eigenstates have been discussed in detail in [15]. And for 
the potential- well-free case (Vq = 0), the eigenfunction of the Hamiltonian Hp = p^/2m — xF 
with (continuous) energy E is given by [TB] 

in terms of the Airy function A\{z). 

The method we employ here to obtain the propagator of the system I-Lf is to apply the 
relevant Lippmann-Schwinger integral equation [HI [IT] 

/CF(x,t|x',0) = iri.(x,t|x',0) + ^ /" rfrKj.(x,t|0,r) ■ /Ci.(0,r|x',0) , (4) 

where the propagator /Cp'(x, 0) := (x| exp(— it As well-known, the propagator 
has the physical meaning of a (complex-valued) transition-probability amplitude to get from 
the point (a;', 0) to the point (x, t) [H]. Here Kp{x, t\x', r) represents the propagator relevant 
to the partial Hamiltonian Hp only, explicitly given by [7j 

Kp{x, t\x', t) = Ko{x - Xc{t), t\x' - x,(r), r} ■ e^yms.{t)-s.{r)} . ^(im{xMt)-.'Mr)} ^ 
in which the propagator 

= \l2mh7t-r) ^"P{^2(r^("-^')'} 

for a free particle subjected to Hq = /2m only. Also, Pc(^) = jldTF{T) = Ft is a field- 
induced classical impulse, leading to the corresponding field-induced translation Xc(t) = 
(1/m) J^drpci^r) = Ft^ /2m and the field-induced action Sc{t) = (l/2m) jldTp1{T) = 
F'^t^/6m. Eq. (|4| immediately gives rise to 

ipp{x,t) = (j)p{x,t) + I dr Kp{x,t\0,T) ■ ^p{0,t) , (7) 

^ Jo 
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where the homogeneous solution 

(ppixjt) = I dx' KF{x,t\x' ,0) %jjo{x') (8) 



obviously represents free motion subjected to an external field F only, while the second 
term on the right-hand side, by construction, gives the infiuence of the residual zero-range 
potential. As an example, for the initial bound state 'il>o{x) = ipbix), Eq. ([s]) reduces to a 
closed expression [7] 

M^,t) = Vb e(V''){-PeW-5cW} .[m(x- x,{ty, -iB- -\ + m( x,{t) - X- -iB- - 

I \ m J \ m 

'(9) 

in terms of the Moshinsky function [19] 



M{x-k-t) = 1 e^(^'^-'='*/2) erfc (^^^^ (10) 

where erfc(2) is the complementary error function |20] . 

It is also instructive to point out that we can employee, as an alternative to the interaction 
Hamiltonian x F in ([T]) given in the scalar-potential gauge, its counterpart pA'm. the vector- 
potential gauge, and the system of interest is accordingly given by 

Ua = ^{P + Pc{t)Y - V,5{x), (11) 

where the vector potential A corresponds to the field- induced impulse Pc{t)- Then it can 
easily be verified that 

KF{x,t\x',T) = exp ^^{xpcit) - x'pc{t)}^ ■ KA{x,t\x',T) (12a) 



iPf{.x, t) = exp\-x- pcit) \ ■ ^pA{x, t) . (12b) 



Substituting (12a) and (12b) into (|4]) and ([T]) respectively, we can straightforwardly obtain 
the equivalent results in the vector-potential gauge to what below follows in the scalar- 
potential gauge (cf. for a detailed discussion oi x F versus pA gauge problem, see, e.g., 
[21]). Besides, our formalism can apply to the system of a delta-potential barrier under an 
electrostatic field 

= 1^ + V^5{x) - xF{t) (13) 

as well, simply by replacing Vq by — Vq in Eq. (|4]) or ([T]) and then going forward straight- 
forwardly. In fact, the eigenstates and eigenvalues of the Hamiltonian I-I'q = Vq 6{x) are 



identical to those of = ~Vo^{x), except for the fact that does not have the bound 
state ipb{x) as its eigenstate. 

To develop the forthcoming discussions in a simplified fashion, let us rescale coordinate 
and time as the dimensionless quantities, x = Bx and t = (hB'^/m) t, respectively [7j, which 
is basically equivalent to setting h = m = Vq = 1. Then the Schrodinger equation for the 
Hamiltonian Hf easily reduces to 

where the relative field strength / = {m/h^B^) F, and Eq. Q accordingly appears as 

ICf{x,i\x',0) = Kf{x,i\S:',0) +i [ df Kf{S:,i\0,f) ■ JCf{0,f\x' ,0) (15) 

Jo 

where the rescaled propagator 



Kfix, i\x', f) = J exp I —J— (x - x'f + + x') (i -r)-'-l-{i- rf 



2Tii{t-T) \ 2{t-fy ' 2' ' 24 

(16) 

which is identical to Kf{x,t — f\x',0) =: Kf{x,x']t — f). For the sake of convenience we 
replace the notation (t, x) by (t, x) for what follows. 

In comparison, we also review briefly the exact results of the field-free case (/ = 0). We 
first set x = in (15), which then enables /Co(0, t\x' , 0) on the left-hand side to immediately 
substitute for /Co(0, t\x', 0) on the right-hand side. Making iterations of the substitution, we 
can finally arrive at the expression [22j 

^ oo 

/Co(x, t\x', 0) = Ko{x, t\x', 0) + - ^(-2;22)""' i"-'erfc(2i) , (17) 



n=l 



where = (|x| -|- \x'\)/\/2ii and Z2 = t/i\/2ii. Here the index n represents a number of the 
substitutions and so that of interactions with the (5-potential well, and i"erfc(z) the repeated 
erfc integral [cf. i°erfc(z) = erfc(z)]. Using the identity |20] 



(7 \ n 
— j e"'erfc(z) = (-2)"n!e"'i"erfc(z), (18) 

we can easily obtain 

oo 

^(-2^2)" i"erfc(zi) = exp (221^2 + z^) erfc {z^ + z^) , (19) 



n=0 



which subsequently leads to the closed expression |22| 

/Co(x,t|x',0) = Ko{x,t\x',0) + ^M{\x\ + \x'\;i]t) . (20) 
Let us come back to the field-induced dynamics (/ 7^ 0) in order to explore an explicit 



expression of ICf{x,t\x' ,0) as the counterpart to the result in (20). For doing so, it may be 



useful to apply the Laplace transform C to Eq. (15) with respect to time; we do this job 



first for X = with the help of the convolution theorem, which will finally give rise to [6] 
C{JCf{x, t\x', 0)}(s) =: Gfix, x'; s) = Gf{x, x'; s) + z GfixAs)Gf{Ox',s) ^ ^^i) 



where the Laplace transform C{Kf{x,x';t)}{s) =: Gf{x,x'; s). From Eq. (16) it can easily 
be shown as well that Gf{x,0;s) = Gf{0,x;s) =: Gf{x]s). To obtain an explicit form of 
Gf{x,x';s), we use the integral representation [7] 



Ai(/3 + \a\) Ci(/3 -\a\) = dt^ exp |^ (^^ - /3t - | , (22) 



where a and P are real- valued, and the Airy function Ci{z) = Bi(z) + iAi{z). This allows 
us to have the Green function in the frequency domain [B] 

/ ~ / / 4 \ ^ TT 

Gf{x,x';s)\s^-iu:+o+ := Gf{x,x';uj) = —7 - Ai(/3 + |a|) Ci(/3 - |a|) , (23) 



,1/1; ^ 

where |a| -> (|/|/4)V3 |x - x'\ and (3 -(2 |/|)-2/3{/ . (3; + x') + 2uj}. Therefore, a 



closed expression of the Green function Gfix, x'; u), identical to Gf{x, x'; s) in Eq. (21 ) with 



s — )■ —iuj + 0+, immediately appears in terms of the Airy functions in (23). However, it is 
highly non-trivial to directly carry out the inverse Fourier transform of the closed expression 
of Qf{x,x']u) in the frequency domain, even with the stationary phase approximation, in 
order to obtain an explicit expression of }Cf{x,t\x' ,0) in the time domain [9]. 

For a later purpose, we introduce here the time-dependent ionization probability for the 
initial bound state 'ipbix), which is defined as 

Pfit) := l-\A^^.{t)\\ (24) 

Here the bound state amplitude 

POO 

A^fit) = / dxMx)^fix,t) = Afit) + Mt)^ (25) 
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where, from Eq. ([T]), 

/oo 
dxipbix) (l)f{x,t) (26a) 
-oo 

/oo 
dxipbix) / dTKf{x,t\0,T)ip{0,T). (26b) 
oo Jo 



Substituting ^ into (26a), we can easily obtain a closed form [7] 



4 _ ,M[if;-i;t) M(-if;-i;t 



which reduces to unity at t = 0, as required. Then it has been shown jEl [71 19] that in the 
weak-field limit |/| ^ 1, the exponential decay law \A^j. (t) p oc e~^f^ is a good approximation 
on the average. This approximation can also be simulated by the ansatz in "smooth" form 

^/(0,r) oc e-'^\ (28) 

where the complex-valued energy E = Ej — ^Tf with Ef = E^ + Af G M and the decay rate 

_ 2 

G M [U E]. Here, A/ is the level shift. The semiclassical value Tj^wkb = e ^i^i is a good 
approximation of for |/| < 1, and Af^wKB = ~|/^ is in excellent agreement to A j up to 
I/I < 0.1. However, this exponential decay approximation, by construction, cannot account 
for the ripples observed in the exact time evolution of the ionization probability resulting 
from the ionization-recombination process in the bound-continuum transition, explicitly 
demonstrated in [7\ from the numerical treatment of the Lippmann-Schwinger equation. 

III. DERIVATION OF AN EXPLICIT EXPRESSION OF PROPAGATOR 

To derive an explicit expression of the propagator, we make the same iterations to Eq. 



(15) as those applied above to the field-free case leading to (17), which then reveals that 

/C/(x,t|x',0) = Kf{x,x';t) + i dr Kf{x,0;t - t) Kf{0, x';t) + 

/o dr Kf{x, 0;t-r) dr' Kf{0, 0; r - r') Kf{0, x'; r') + dr Kf{x, 0; t - r) x 
/; dr' Kf{0, 0; r - r') dr" Kf{0, 0; r' - r") Kf{0, x'; r") + ■ ■ ■ . (29) 

To simplify the notation, let Kf{x, 0; t) = Kf{0, x; t) =: Kf{x; t) and then -^'/(O; t) =: Kf{t) 



from now on [cf. (16)]. For a later purpose, we rewrite Eq. (29) as 



]Cf{x,t\x',0) = Kf{x,x']t) + i [ dTKf{x;t-T) [ dr' Kf{x';T - r') ■ A{t') , (30) 

Jo Jo 



where the partial integrand, being a function of time only, 

A(r') := 6{t') + iKf{T') + f dr" Kf{r' - r") Kf{r") + 

Jo 



r dr" Kj{r' - r") f dr'" Kf{r" - r'") Kj{t"') + ■■■ . (31) 
Jo Jo 



Eq. (30) can easily be verified with the help of the Laplace transform in (21), but without 



considering Eq. (23); the Laplace transform is straightforwardly expanded as 



Gfix, x'; s) = Gf{x, x'- s) + I Gf{x; s) Gf{x'; s) ^ {z Gf{0; s)}" . (32) 



n=0 



And then we simply apply the inverse Laplace transform to this expression term by term. 



which readily recovers the result in (30). From this, it also follows that 



^f{x,t) = (j)f{x,t) +i / dTKf{x;t-r) / dr' (j)f{0,T - t') ■ A{t') . (33) 

Jo Jo 

As seen, the finding of an explicit expression of the quantity A(r') is a critical step to deriving 
the wavefunction ^f(x,t) for any initial state i^o^x) [cf. Eq. ([T])]. 



Comments deserve here. First, the substitution of x = x' = into Eq. (30) and then the 



comparison with (31) gives rise to a noteworthy relation 

/C/(0, t|0, 0) = i 5{t) - i A{t) , (34) 



which will be useful later. Secondly, the quantity Gf{0;s) in (32) designates the Green 
function for a closed path with the initial and final point being x = 0, and the summa- 
tion index n represents how many times the closed path interacts with the potential well 
located at x = 0. Then the quantity A(r'), exactly corresponding to this summation in the 
time domain, is accordingly responsible for the ionization-recombination resulting from the 
scattering from the zero-range potential at x = 0, particularly for tpo{x) = ipb{x), which 
was qualitatively explained in [TJ [13] as follows: This initial bound state, explicitly given by 



'■Pbip) = a/2 /tt / (j9^ + 1) in the momentum representation, has a symmetric momentum dis- 
tribution around p = for the particles (in an ensemble represented by the bound state). By 
applying the field, this symmetry breaks down in such a way that the motion of the particles 
in one direction is accelerated and so they will easily leave the potential well, simply toward 
the (continuous) unbound states. On the other hand, the motion in the other direction 
gets slowed down until the particles stop, and then they reverse their direction of motion so 



that the particles again approach the potential well at x = where they cause a transient 
maximum in the bound state probability. This process is repeated until all particles will 
completely leave the potential well. Below we will systematically derive the quantity A(r'), 
and so /C/(0, t|0, 0), in closed form. 

Let first Ao(t') := 6{t') and Ai(r') := iKf{T') in (31). We consider the next term 
A2(t') = — I2(t'), where the integral 



hit) 



drKf{t-r)Kf{r) 



dr 



27ri Jq y/t -t y/r 



exp 



24 



{{t-rr + r'}) (35) 



[cf. (16)]. Let r = ty. Applying to Eq. (35) the Taylor series of the exponential function 



and then the identity of beta function 



\yy-Hl-yr' = ^4^ = Bia,p) 



r(a + /3) 



where Re(a),Re(/3) > 0, it easily appears that 
hit) 



27ii ^ k\l\\ 2Ai 

k,l=0 



(36) 



(37) 



r(3fc + 3/ + i) 

With the help of the Gamma function identities, T{3z) = {2n)-^ 33.-1/2 y(^z) T{z + 1/3) T{z + 



2/3) and 1(1/6)1(1/2)1(5/6) = 2tt^/^ and 1(1/3)1(2/3) = 2n/3'^/^ l20j, Eq. ^ reduces 
to a closed expression 

. ri 1 51 . n 1 51 

• L6' 2' eJ ' L6' 2' eJ 

LS' 3' J ■ ' 

with a = f'^t^/24i in terms of the generalized multiple hypergeometric function 



hit) 



— X F ^- 3; 3 
2i 3:-;- 



a, a 



(38) 



pPo-Ph- ■ 




00 /X ^ ( \ 



k,,-,k„=o (bo)fci+...+fc„ fJi {hj)k, f^j 



(39) 



where = (a^i, ■ ■ ■ ,ajpj) and hj = {bji, 



with j = 0, 1, ■ ■ ■ ,n are vectors with 



dimensions pj and qj, respectively. And (aj)^ := H/'iil'^iOfe i.^j)k '■= Y[lLi{bji)k, where 



the Pochhammer symbol (A)^ = r(A + A;)/r(A). The multiple series in (39) absolutely 
converges if 1 + go + Qj' — Po — Pj' ^ for all j' = 1, 2, ■ ■ ■ , n. In fact, Eq. (38) fulfills this 
condition as 1 + 3 + - - 3 = 1 > for j' = 1, 2. 



Subsequently we consider the next term A^i^r') = —il^i^r'), where the integral 



hit) 



drKfit 



r 



dr' Kfir-r') Kfir'). 



(40) 
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Similarly to l2{t), we can straightforwardly obtain 

^ oo / ^2 \ k+l+n 



/c,/ ,Tl — 

(41) 

Let r' = ry' and r = ty. We then use the identity [24 



' du rft; .."-^ t;^-^ (1 - - vV-^ = , (42) 

' ^ ' r(a + /3 + 7) ' ^ ' 



Jo 



where Re(a;), Re(/3), Re(7) > 0. Eq. (42) is easily shown to be identical to JJ" rfw (1 — 

t/;)'^""'^ {w — u)'^~'^ u""-*^, where w = 1 — v and J^^ rfw dw = dw du. Consequently it 
appears that 

^ r^l^ V ^ r(3fc + ^)r(3/ + ^)r(3n+^) 

^'^^ \(27r03/ ^Z^^A;!/!n! V242y' r (3A; + 3/ + 3n + |) ' ^ ' 



With the aid of the technique already applied to (37), this reduces to a closed expression 

T (f\ - f ^ \ ' F 3; 3; 3 ( ~ ' ts' 2' e] ' 2' e] ' 2' e] ' ^ ^ ^ | (aa\ 

\ L2' 6' eJ ■ ' ' ' / 

with a = pt^ /24:i. Along the same line, we can subsequently obtain closed expressions of 



l4(t) and l5(t), ■ ■ ■ , respectively, where ^^{t') = l^ir') and A5(r') = i I5(t'), ■ ■ ■ for Eq. (31 ). 

Now let us generalize the above scenario to the quantity A„(r') = 2"I„(r'), where the 
integral 

I„(t) := / dTiKf{t-Ti) dT2Kf{Ti-T2)--- dTn-lKf{Tn-2- Tn-l) Kf{Tn^i) . 

Jo Jo Jo 

(45) 

To do so, we apply the same technique with the aid of the identity |24] 

pl—ui rl-ui Up-i 

dui / du2 ■■■ dUp ■ ■ ■ UpP'^ (1 - Ml Mp)^"^ 



^ r(aO---r(a,)r(/3) 

r(ai + ■ ■ ■ + Op + ^) ' ^ ^ 

where Re(/3), Re(aj) > for j = 1, 2, ■ ■ ■ , p. This allows us to finally arrive at the expression 

1 f^y (/^t3/24.)^^-^-+^" r(3fc, + i)---r(3A;„ + i) 

"^^ t\2mj ^^^^^^^ h\---u r{3h + ... + 3K + ^) ' 
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which subsequently reduces to the closed expression 



In(t) 



r(f) 

: 3;---: 3 



X 



(48) 



. n 1 51 
• Le' 2' eJ 



n 1 51 . 
Le' 2' eJ ' 



fn ra+2 n+4 1 . _ 
Lg' 6 ' 6 J • (1) 



(n) 



As a result, we find the exact expression 

oo 



(49) 



n=2 



in terms of the (well-defined) multiple hypergeometric functions. From (34), it easily follows 
as well that /C/(0,t|0,0) = Kf{t) + E^=2 ^""^ ^"(i)- By substituting Eqs. (|l6| and (g 



into (30), we can now obtain an explicit expression of the propagator /Cj(a;, 0) in terms 
of the single double integral dr dr' {■ ■ ■ A(r')}; this remaining double integral cannot 
straightforwardly be evaluated indeed in terms of the multiple hypergeometric functions in 



(39). However, it is still desirable to express Eq. (49) (in good approximation) in terms 
of the generalized hypergeometric functions pFg{z) (with a single argument z) being 
much more transparent to physically interpret and much more accessible in their numerical 



evaluation than the multiple hypergeometric functions [cf. (83) and (84) 



In order to discuss the mathematical complexity of Eq. (49), it is also instructive to 
formally express the quantity A(r') in terms of a single integral over the ri-sphere solid angle 
Qn (e-g-; a 1-sphere corresponding to a circle), in which the infinitesimal element dfln '■= 
(sin 0n)"~^ ■ ■ ■ sin 02 d(f)n ■ ■ ■ d(f)2 d(j)i, and the integral over the entire region is explicitly given 
by § dVLn = Jq d(f)n {sin (f)n)"'~^ ■■ ■ /J^ (i02 sin 02 Jq^ d(l)i [27j. To this end, we first consider 
the integral l2(t) in (35). Let r = t(sin^)^. This immediately allows us to have 



hit) = - 
m 

which reduces to a closed form 

hit) 



dO exp 



24 



{(cos^)^ + (sin^)^} 



(50) 



2i 



exp 



192 



PI 
64 



in terms of the Bessel function J^iz) [cf. (38)]. Here we used (cos 6*)^ + (sin( 
(3/8) • (cos 4^) and then J^iz) = (l/vr) d(f) cos(;zcos0) [20]. 



(51) 

5/8 + 
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The integral l3(t) is next under consideration. We substitute Eqs. (16) and (51) into 



(40) and then carry out the same change of integral variable as that used for (50). Then it 



appears that 



hit) 



t 



2ttP 



dO sin 6 ■ exp 



{8 (cos^)^ + 5 {sinef}] Jo 



192 



64 



(52) 



[cf. (65) for exact evaluation of this integral]. We next pay an explicit attention to 



Uit) 



dn Kf{t - n) / dr2 Kfin - r^) / dn Kf{r2 - rg) Kfir^) (53) 



[cf. Eq. (45) with n = 4]. Using dr JJ^ dr' = dr' J^, dr and then introducing u := r — r', 
the first two-integral part of ( 53 ) is easily transformed into dr' J^^ ^ duKf{t—T'—u)Kf{u). 
This, with Eqs. (35) and (51), then allows us to have 



dr exp 



192 



fit 



64 



64 



(54) 



Applying again the very technique already used for l2(t) and l3(t), we can finally obtain 



dd (sin 2^) exp 



2+3 



/^i(:^(cosi 

6r~ 



Jn 



192 
64 



■{(cos^)^ + (sin^)^} X 



(55) 



[cf. (76)]. Similarly, we can also find that 



hit) 



S,) ^ A (sin 9)^ (sin 2^) Jo ( m^iE^lf^ 



J I ft^jsine)'^ (siny)*^ 



X 



2+3 



exp<' (8 (cos^)^ + 5(sin0)6{(cos(/^)'' + (sin(/^)'^}) 



(56) 



192 

where dQ2 ■ ■ ■ = Jq^"^ f^^"^ sin 9d9dip ■ ■ ■ represents an integral over the surface of a unit 
2-sphere (with radius r = 1) covering the first octant only. Along the same line, we continue 
to be able to obtain the corresponding expressions of l^it) and l7(t), ■ ■ ■ , respectively. In 
doing this job, it can easily be induced that 



hm+iit) 

hm+2it) 



2t\2 

Tci J 



sm (pm ■ exp 



t 



d(j)m (sin 2(j)m) exp 
hm (it(sin0m)^) , 



—if r (cos( 

5ipt^ (cos (prr. 

192 



Jn 



l2m (t (sin0m)^) (57a) 

(57b) 



64 
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where m = 1, 2, • • • . 

Now let us generalize the above scenario to the integral ln{t). Applying the iterations to 



Eqs. (51), (57a) and (57b) allows us to finally obtain the following irreducible expressions 



in terms of an integral over the higher-dimensional solid angle as 



2nit / V i 



ifH^ (COS0^)^ 

dilm exp <j ^ > X 



n 

k=l 



[SVClCpkt {cOS(t)k~l) exp ( [COS (pm-k) 



192 



n 



Ism (pi) 



l=m—k+l 



Jo I {cos (j)m-kf 



n 



ism (pi 



l=m—k+l 




(58a) 



and 



m+l 



I 



2m+2 



it) 



2i\i) Jl (sin 0^+1)™+! 




l=m-k+2 / \ ^ ^ ^ ' l=m-k+2 

Here the angle (pQ := 0, and j\^dQm--- represents an integral over the first section only, 
corresponding to < ,(pm ^ 71/2. Now we easily see that with the field strength 

/ — > 00, the integrals ln{t) — )■ for n > 2, due to the fact that Jo{z) — )■ with z — )■ 00. As a 



result, the quantity A(r') in (49) can be rewritten in terms of the integrals over the higher- 



dimensional solid angle. This result can also be interpreted in such a way that Eq. (31) 
for A(r'), given in form of the time-ordered integrals dr" dr'" ■ ■ ■ and so non-trivial 
to directly evaluate, is transformed, for an arbitrary time r', into the integrals of some 
geometric pattern over the (time-independent) surface of a unit n-sphere, so being more 



accessible in numerical evaluation; in fact, we see from (58a) and (58b) that the quantity 



I„(t) is given by an irreducible [{n — l)/2] -dimensional integral, rather than an {n — 1)- 



dimensional integral in (31), where the symbol [y] is the greatest integer less than or equal 
to y. Due to this irreducible high dimensionality, it is, apparently, highly non-trivial, though, 
to exactly evaluate the (complicated) integral In(t) for n being large enough in terms of the 
hypergeometric functions pFq{z) simply with a single argument z. Below we will accordingly 
explore an approximation scheme in which the quantity A(r') can be expressed in terms of 
the functions pFq{z) in reasonably simple form. 
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IV. APPROXIMATION OF PROPAGATOR BASED ON THE PARTIAL WAVE 
EXPANSION 



We already have the closed expression of the integral l2{t) in terms of the Bessel function 



in (51 ). We next intend to evaluate the integral Isit) in (52) in terms of some hypergeometric 
function. Substituting first the identity |20] 



JJz) 



[z/2Ye-'^ ( 1 



(59) 



with z/ = into (|52|), we can straightforwardly obtain 

t 



where 



l3l(t) 



exp 



192 



n=0 



32 



de (sin0)^"+^ exp ( -^4^ cos 4^ 



(60) 



(61) 



'0 V 64 

Here the confluent hypergeometric function iFi{a; b; z) = {r(6)/r(a)} Yl'^=oi^(^~^^)/^(^~^ 
n)} z^/n\. 

Now we focus on explicit evaluation of the integral l3i(t). To do so, we flrst substitute 



into (61) the expansion formula for a plane wave |20] 

oo 

exp {ia COS (j)) = Ji{a) cos(/</)) 



(62) 



l=—oo 



with — )■ 4^ and a where the Bessel functions Ji{z). As seen, each / of the 

harmonic partial waves is then time-independent. Subsequently, with the help of another 
sum rule cos(?t,z) = Ylt=o iT) (cos 2;)'^ (sin 2;)""'^' cos{(n — A;) 7r/2} with n = Al followed by the 
beta function B{w,z) = 2 JJ"^^ (sin0)^"'"^ 
integral representation |i28j 



cos ( 



\2z-l 



[cf. (36)], we can easily obtain the 



de(smey(cosAW) = . — ^\ ^ 

^ ^ ^ ^ 2P+1 r (1 + 2/ + 2) r (1 - 2/ + 1) 



(63) 



in which p ^ 6n + 1. In doing so, we also used r{z)r{l — z) = 7rcsc(7r2;) and T{2z) 



(27r)-i/2 22^-1/2 Y(^z) V{z + 1/2) [20j. Eqs. (g and (g then allow us to have an evaluation 
of the integral l3i(t) and so 



h{t) = {^Y 



r(i+«).(6n+i)! ( ift^ Y Sr _ 



n=0 



Z=— 00 



f+3n+2«jr(f+3n-2/j 



(64) 
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This subsequently reduces to a compact expression 



hit) 



2nr 



1/2 



exp 



192 



l=—oo 



- 16/2] 



Jj 



PI 

64 



(6) (^Pt' 



32 



(65) 



in terms of the generahzed hypergeometric functions 



Z) := 6^6 



1 1 1 2 5 7' 



3' 2' 2' 3' 6' 6 



2/1 2/5 

3"' 2 ^ 3"' 6 



2/5 2/7 
3"'6 ^ 3"'6 



2/ 7 2/ 

3"'6 ^ 3" 



which satisfies 

oo 

E 



r (I +n) ■ (6n+ 1)! 



n=0 

Here we also used JA—z 



r(l + n) r (I + 2/ + 3n) T (| - 2/ + 3n) 



n! 



4-7r-V2 
1 - 16/2 



f/^) (2«.) 



(66) 



(67) 



'-YMz). 



Comments deserve here. First, Eq. (65) can be rewritten as 



l3(t) 



t 

2^ 



exp 



Et 



2i- 



1=1 



16/2 




5 7' 
6' 6 



^pf 
32 



+ 



(68) 



[cf. (44) and (59)]. Now it is explicitly shown that the quantity l2{t) in (51) is "modulated" 



by 2^2(' ■ ■ ) yielding in (68) the partial wave / = as the leading term of l3(t), surrounded by 
the additional partial waves / ^ "modulated" by f/^^(- ■ ■ )• Secondly, we may accordingly 
take the leading term only as a satisfactory approximation of Isit); in the weak- field regime 
(/ ^ 1) leading to Ji{b) —i- for / 7^ with b = f'^t^/6A, we easily see the validity of this 
approximation. It also applies sufficiently to the strong-field regime (/ ^ 1), in which due 
to the asymptotic behavior Ji{b) ^ {2 / (nb)}'^/'^ cos{b - Iti/2 - tx/A) [20J, the magnitude of 
all terms with / 7^ cannot be non-negligible enough. In fact, we have |l3(t)| <^ 1 anyway in 



this regime, as already pointed out after Eq. (58b). Also, the strong- field regime corresponds 
to the semiclassical limit /i — i- since the argument 6, expressed in a dimensionless unit, 
exactly corresponds to F'^t^/{64hm) in the actual physical unit [cf. (14)]. Therefore the 
approximation with / = alone may also be considered an effective semiclassical treatment 
of hit). 

Next we explore a closed expression of l4(t) in approximation. The substitution of the 



leading term of l3(t) into Eq. (53) allows us to straightforwardly obtain 



hit) 



1 ^ T{^ + n) r(2 + 6r2) 

'^h{ni+n)r{i+3n)y 



^_P 
211 



T{\ + k) [if 



- {1(1 + ^)1^2^ 



xl4i(t), (69) 
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where 



l4l(t) 



t ^3{n+k)+l/2 / 

dr , exp I 

V 24 



2 ^3(n+fc)+l 



exp 



192 



dO (sin^)^("+'=)+2 exp 







64 



cos 4^ . (70) 



Here we applied the same technique as that used for (60 )-(61 ). Subsequently we again apply 



to Eq. (70) both the partial wave expansion in (62), followed by the selection of / = alone 



and Eq. (63) with p — Q{n + A;) + 2. This immediately gives rise to 

3(n+A:)+l 



Ulit) 



exp 



^ipt^\ J fft^\ T{6{n + k) + 3} 
Jo 



192 



64 J (r{3(n + A;) + 2})2 



(71) 



Substituting this into (69), we can next obtain straightforwardly 



192 



^0 ( ^ ) X l42(t) 



(72) 



where 



_ ^ r(l+n) r(2 + 6n) f^fH■^Y^T{l + k)^{3 + 6{n + k)} frfH^^' 



n=0 



{r (I + 3n) -niy 



k=0 



(r{2 + 3(n + A;)}-A;!)^ 



(73) 



Let n + k = r. Then, this can easily be rewritten as 



r(3 + 6r) f^fH^y^ T{l+n)T{2 + 6n)T{l + r-n) fiy 



26 



.=o--'-- ^ n=o {r(| + 3n) r(l + r-n)-n!} 

Here the second summation over the index n precisely simplifies to 4 ■ F + r) /{^/tt (r!)^} • 

iih i' ~'"] ' [i' 6' ^ ~ '"l ' ~''') ^^'-^ '^^^ approximated with satisfactory precision 
to its leading term n = alone, due to the fact that the summand (1/2^)" — > for all n > 1 
(cf. Fig. [T]). From this, Eq. (74) reduces to 



l42(t) 



_4_ ^ F(^+r) F(3 + 6r) ft ft 
{F(2 + 3r)-r!}' 



11 



84F4 



115 7' 
2' 2' 6' 6 



3' ' '3 



32 



which immediately gives rise to the expression 
t 



hit) 



exp 



-'0 1 I 4^4 



192 



64 



115 7 
2' 2' 6' 6 



-,1,1, - 
3' ' '3 



32 



(75) 



(76) 



in terms of the generalized hypergeometric function. As demonstrated in (69)- (75), without 



this leading-term approximation the quantity l^it) should be expressed in terms of a lengthy 



17 



multiple sum, which is not transparent to physically interpret and not desirable for numerical 
evaluation. 

Next the quantity l5{t) is under consideration. We exactly follow the technique leading to 



Eq. (76) for l4(t); we first substitute this previous result with (75) into (45) with n = 5 and 



then carry out the same approximation, followed by applying Eq. (63) with p — 6(n + A;) + 3 
From this, we can arrive at the rather lengthy expression 

^3/2 



1 



64 V 2m 



exp 



192 



64 



(77) 



where 



^^^^^^ __^ r{l + n) r(3 + 6n) / ^fHy ^ TC^ + k)^{A + 6{n + k)} ^ ^fHy ^^^^ 



n=0 



{T{2 + 3n)-n\y 



k=0 



{T {I + 3{n + k)} ■ k\y 



cf. (73)]. We now apply the technique used for (74) and (75), which finally gives rise to 

^3/2 



hit) 

where 



1 

2^ 



exp 



-'o ^rr 3^3 



192 



64 



32 



E 



r(l + r) r(4 + 6r) fifH 



t {T{l + 3r).r\y 



)11 



3^3 



1 2 4 
2' 3' 3 



1 2 4 

2' 3' 3 



5 7 3 
6' 6' 2 



5 7 3 
6' 6' 2 



32 



32 



(79) 



(80) 



[cf. (75)]. Along the same line, we continue to be able to obtain the next expressions in 



approximation as 



hit) 



it' 



16 V 192i 



-'0 I -^ry- I 4fM 



64 



15 7 3 
2' 6' 6' 2 



1,1, 



4 5 
3'3 



32 



(81) 



hiit) 
Ii2(t) 



{i2m)-H^}-: 
3 ■ 5 ■ 7 ■ 9i 

exp 



exp 



26-5! 



192i 
5^/2^3 



64 



4F4 



2' 3 3 



192 



1 11 13 5 
2' Y' ~6' 2 



11 13 5 
IT' IT' 2 



32 



i.2.;« 

' '3'3 



32 



By induction it easily follows that 
hit) ^ hit)-^l^] 4F, 



r (f ) V2^ 



1 n — 1 n + 1 n + 3 
2' 6 ■ 



6 



6 



n n + 2 n + A 



"'6' 6 



6 



^Pt' 



32 



(82) 



where n = 3,4, ■ ■ ■ [cf. Eqs. (58a) and (58b)]. Here the expression with n = 3 is obviously 



meant as the leading term of Eq. (68) only. It is also instructive to note that letting / — )■ 0, 
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Eq. ( [82) ) directly recovers the exact expression of its field-free counterpart, explicitly given 
by - l)/{(2i)'^/2 . T{n/2)}, obtained from substitution of ^ with / = into Q. 
Now we are ready to have a closed expression of the quantity A(r') in approximation, 



accounting for the ionization-recombination, as discussed after Eqs. (|28|) and (|33|). The 



substitution of (82) into (49) explicitly reveals that 



Aa(r') ^ 6{r') + iKfir') + - exp 



192 



64 



4F4 



1 n — 1 n + 1 n + 3 
2' ■ ■ 



6 



6 



6 



n n + 2 n + 4 if t 



n=3 
2^/3 



r(i) 



IT 

~2 



X 



6' 



6 



6 



32 



(83) 



This expression is useful especially in the weak-field limit, as stated after Eq. (82). To see 



transparently the strong-field behavior of Aa(r'), it may also be desirable to rewrite (83) as 
follows; we first decompose the sum index n as 6A; + 3, 6k + A,6k + 5, 6k + 6,6k + 7, and 
6k + 8, where k = 0,1, ■ ■ ■ . Next we plug the expansion ■■];[■■■]; z) = Yl'^oi' ' ' ) -^V^' 
into (83). With the aid of the sum identity X^fclo Xl/^o 5'(^' — Xl^^o X]I=o 5'(^' ^ ~ ^) 
where r = k + I, we can unify and then simplify the two r'-dependencies, explicitly given 



by (■ ■ ■ r')"/2"^ and 4F4{ 
Aa(r') as 



in the infinite sum over n. This finally allows us to obtain 



h(T') i 1 + 



f 



I 5 



j=o \ 2 J 



IT 

~2 



E 

fc=0 



r(l + | + 3A;) r(A;+i) {r(l - A;)}^ 



X 



5^5 



2+j 4 + J 6 + j 1 



6 



6 



6 ' 2 



3 + J 5+j 7 + j 



kA-k 



2^/3 



32 



4) 



6 ' 6 ' 6 ' 

with / — !■ 00 [note that the sum over j is just a finite one]. Here the generalized hypergeomet 
ric function ^^{- ■ ■ i/^r'^/32) can be understood as a compact form of Re(- ■ ■) +i Im(- ■ ■ ) 
explicitly given by 

1 - 2A; 3 - 2fc 2 + j 4 + j 6 + j 8 + j 10 + j 12 + j 



9^] 



^° ' 4 ' 4 ' 12 ' 12 ' 12 ' 12 ' 12 ' 12 

l-k l~k 2-k 2-k 3 + j 5 + j 7 + j 9 + j 11 + j 13 + j 



2 ' 2 ' 2 ' 2 ' 12 ' 
/V^ /1-2A;\ /2+j\ /4 + j 



32 
9^10 



6 



6 



12 ' 12 

6+J 
6 



12 



12 



12 



)12 



X 



3-2A; 5 -2k 8 + j 10 + j 12 + j U + j 16 + j 18 + j 
' 4 ' 4 ' 12 ' 12 ' 12 ' 12 ' 12 ' 12 
l-k l-k 2-k 2-k 3 + j 5 + j 7 + j 9 + j 11 + j 13 + j 



12 



12 



12 



12 



12 



12 



4^/6" 



)12 



.(85) 
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And it is worthwhile to comment that although the function ^F^il- ■■];[■■■ A — k,l — k]; ■ ■ ■) 
itself is not well-defined for k = 1,2, ■ ■ ■ , the composite expression ^F^d- ■■];[■■■ , 1 — ^, 1 — 



k]; ■ ■ ■ )/{r(l — fc)}^, given in (84) indeed, is mathematically undisputed. As a result, the 



substitution of the quantity Aa(r') in (83) or (84) into (30) can directly yield an explicit 
expression of the propagator 

ICa{x,t\x,0) ^ Kf{x,x';t) + it^ / dyyKf{x;t{l-y)) / dy Kf{x';ty{l - y')) ■ Aa{tyy') 



(86) 

in terms of the generalized hypergeometric functions and an integral which can straight- 
forwardly be evaluated to any sufficient degree of precision. Here we used r' = ry' and 
r = ty. In fact, it can easily be shown that this remaining double integral is beyond the 
scope of its evaluation in terms of any generalized hypergeometric functions in reasonably 



simple form. From (33), we can also obtain the corresponding wavefunction ijja{x,t) for an 
arbitrary initial state. 

Comments deserve here. The wavefunction il)a{x,t) evolved from the initial bound state 
ipb{x) must, by construction, accommodate the ripples observed in the ionization probability 
Vf{t) in (24). To have this feature, we took the leading term only from each integral In(t): 



In the strong-field limit, the first term (f)f[x,t) on the right-hand side of Eq. (33) for 
il)a{x,t) is dominant to the second term, given by a double integral, which has therefore 



been neglected in the analytical approach in [7J. Now, the result in (86) allows us to study 
systematically the next-order terms to (f)f{x,t) (in the intermediate-field regime), which are 
responsible indeed for the ripples resulting from the scattering from the delta-potential well 
at X = 0. In the weak-field limit, on the other hand, the second term (corresponding to the 
infiuence of the residual zero-range potential) is dominant to the first term subjected to the 
external field only, and so the quantity Aa(r'), and so /Ca(0, t|0, 0), is a highly critical factor 
to determination of the entire time-evolution ipa{x,t). In fact, the bound-state amplitude 
A^j{t) = A^j,{t) + As{t) of the ionization probability in Eqs. (24)-(27) then reduces to a 
compact form as 

As{t)-^t[ dr [ dT'<j)f{0,t-T)<f)f{0,r-r')Aa{r') (87) 
Jo Jo 

[cf. ([9]) for the closed expression of 0/(0, t)]. Finally, we also remark that in actual numerical 
evaluation of (86) we need to introduce a (small) constant c(t) in order to exactly fulfill the 
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normalization condition dx \''pa{x, t)p = 1 in such a way that ipa{x, t) = (pf{x, t) + c(t) x 
it'' Jldy j^dy' {■■■). 

V. CONCLUSION 

In summary, we have investigated the time evolution of a particle subjected to both a uni- 
form electrostatic field and an attractive one-dimensional delta-function potential. We have 
systematically derived the propagator /C/(x, 0) of this system, in which its essential in- 
gredient /Cj(0, )f:|0, 0), accounting for the ionization-recombination in the bound-continuum 
transition, is exactly expressed in terms of the multiple hypergeometric functions. And 
then, as our central finding, we have obtained the ingredient /C/(0,t|0,0) in (appropriate) 
approximation, expressed in terms of the generalized hypergeometric functions being much 
more transparent to physically interpret and much more accessible in their numerical eval- 
uation than the multiple hypergeometric functions. It has been shown that this finding 
provides a much better approximation scheme than the exponential decay approximation, 
where no ionization-recombination dynamics can be treated. In our approach, we have not 
applied the energy-time Fourier nor the equivalent Laplace transform, which has, on the 
other hand, been normally a starting point for analytical study of this system in references, 
e.g., |Sl m El 123] ■ We think that our approach will therefore provide a useful framework for 
analytical study of the time-dependent ionization process in a delta-function potential well 
under an oscillatory electric field as the next task, for which one cannot straightforwardly 
apply the convolution theorem of the Laplace transform to the relevant Lippmann-Schwinger 
integral equation. 

Acknowledgments 

The author would like to thank G.J. lafrate (NC State University), Y. Wang, and L. 
Zhou for stimulating discussions. He also appreciates support from G.J. lafrate during his 



21 



stay with NC State University. 



[1] T. Yamabe, A. Tachibana, and H.J. Silverstone, Phys. Rev. A 16, 877 (1977), and references 
therein. 

[2] M. Razavy, Quantum Theory of Tunneling (World Scientific, Singapore, 2003). 
[3] C. Joachim and J.R Launay, Chem. Phys. 109, 93 (1986). 

[4] M. Jonson, Phys. Rev. B 39, 5924 (1989). 

[5] S. Geltman, J. Phys. B 10, 831 (1977) and J. Phys. B 11, 3323 (1978). 

[6] A. Ludviksson, J. Phys. A 20, 4733 (1987). 

[7] W. Elberfeld and M. Kleber, Z. Phys. B 73, 23 (1988). 

[8] S.M. Susskind, S.C. Cowley and E.J. Valeo, Phys. Rev. A 42, 3090 (1990). 

[9] M. Kleber, Phys. Rep. 236, 331 (1994). 
[10] B.-G. Englert, Lett. Math. Phys. 34, 239 (1995); here, an "analytical expression" of the ion- 
ization probability Vf{t) was derived, which is, however, just a complicated formal expression 
in terms of some differential operator [cf. Eqs. (29) and (35) thereof]. 
[11] A. Rokhlenko and J.L. Lebowitz, J. Math. Phys. 41, 3511 (2000); here, the time-dependent 
wavefunction was expanded in terms of a complete set of eigenf unctions of the field-free Hami- 
tonian Ho = p^/2m — VoS{x) to study the system of a single delta-function potential well 
which is subject to time-dependent variations of the binding strength, but with no external 
field. 

[12] W. Elberfeld and M. Kleber, Ionization and tunneling in a strong electric field, in Physics of 

Strong Fields, edited by W. Greiner (Plenum, New York, 1987) p. 465. 
[13] W.P. Reinhardt, Time- dependent approach to spectra in external fields, in Atomic Excitation 

and Recombination in External Fields, edited by M.H. Nayfeh and C.W. Clark (Gordon and 

Breach, New York, 1985). 
[14] D.J. Griffiths, Introduction to Quantum Mechanics, 2nd ed. (Pearson Prentice Hall, New Jersey, 

2005). 

[15] W.C. Damert, Am. J. Phys. 43(6), 531 (1975). 

[16] O. Vallee and M. Soares, Airy Functions and Applications to Physics, 2nd ed. (World Scientific, 
Singapore, 2010). 

22 



[17] S. Weinberg, The Quantum Theory of Fields (Cambridge University Press, 1995). 

[18] C. Grosche C and F. Steiner, Handbook of Feynman path integrals (Springer, New York, 1998). 

[19] M. Moshinsky, Phys. Rev. 88, 625 (1952). 

[20] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, 

and Mathematical Tables (Dover, New York, 1974). 
[21] R.R. Schlicher, W. Becker, J. Bergou, and M.O. Scully, Interaction Hamiltonian in Quantum 

Optics or p A vs. Er Revisited, in Quantum Electrodynamics and Quantum Optics, edited by 

A.O. Barut (Plenum, New York, 1984) p. 405. 
[22] D. Bauch, Nuovo Cimento B 85, 118 (1985). 
[23] W. Elberfeld and M. Kleber, Am. J. Phys. 56(2), 154 (1988). 

[24] H. Exton, Multiple Hypergeometric Functions and Applications (Halsted/Wiley, New York, 
1976). 

[25] H.M. Srivastava and P.W. Karlsson, Multiple Gaussian Hypergeometric Series (Halsted/Wiley, 

New York, 1985). 
[26] H.M. Srivastava, J. Phys. A 20, 847 (1987). 

[27] D.M.Y. Sommerville, An introduction to the geometry of n dimensions (Dover, New York, 
1958). 



[28] Eqs. ( |63[ ), (65) and (67) as well as all following expressions in terms of the hypergeometric 
functions have been verified with the help of Maple programs. 



Fig. |T| (Color online) y = AV{\+r) (r!)^} 4F3([|, |, -r, -rj; % \, \ - r]; -1) versus 

index r = 0,1, 2,--- with an interpolated line (red dash). This precisely represents the 



second summation over the index n in Eq. (74). In comparison, its approximation ya = 
r(|)r(2)r(^ +r)/{r(|)r(l +r)}^ with an interpolated line (blue dash), where y > ya- 
This is simply the leading term n = of the summation. As seen, y y^ m. a, satisfactory 
manner. 
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